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1. Introduction 



It has been know for some time that electrons on a plane in the presence of a strong magnetic 
field perpendicular to the plane are a realization of a noncommutative geometry. One such 
realization was proposed by Susskind[|l| who describes the quantum Hall fluid by a noncom- 
mutative Chern-Simons theory. Susskind's theory is an effective theory which is lowest (first) 
order in derivatives. However it describes only the quantum state of a given conductance 
plateau. It does not seem to contain sufficient dynamics to describe transitions between lev- 
els, nor the end transition of any given sample to the Hall insulator. We imagine that such 
a theory would contain at least two different kinds of excitations, for example plane waves 
and solitons. In one region of parameter space the plane waves would be the light and hence 
dominant mode while in another it would be the solitons. As a function of a parameter, 
say the magnetic field, if the two excitations became degenerate, then one would expect a 
transition in the behaviour of the system at that point. 

A pure Chern-Simons theory does not contain excitations. A possible second order gen- 
eralization corresponds to adding a Maxwell term to the action. Such a modified theory 
contains a richer dynamical content which may allow for a description of the transitions. A 
version of noncommutative Maxwell-Chern-Simons theory has been shown to contain vortex- 
like solitons ||, however plane wave solutions have not been found there. In that article the 
original Chern-Simons action described by PolychronakosjiJ] was used. This action yields 
homogeneous equations of motion and it is quite straightforward to find vortex-like solitons 
once the vacuum configuration has been found. This version of Chern-Simons theory does not 
admit a smooth commutative limit. In a second version of noncommutative Maxwell-Chern- 
Simons theory, indeed the version studied by Susskind appended by the Maxwell term, it is 
possible to find plane wave solutions E|. However solitons do not seem to exist in this theory. 
It seems that the vacuum of the theory is unstable and solitons are possibly incompatible 
with an infinite plane. However in the analysis there|2| an intriguing possibility suggested 
that they might be present in finite quantum Hall droplets. In this paper we study exactly 
this theory of finite quantum Hall droplets and confirm the existence of soliton solutions. We 
find a rich structure of soliton solutions and we compute their properties and energies. 

2. Action and equations of motion 

Following Susskind's Q notation we describe the quantum Hall fluid with two fields: 

x i = y i + 6e ij Aj Vi e {1, 2} (2.1) 

Ai parametrizes the deviations from the equilibrium situation, x % = y l , which is a static 
solution {y l = = x % is the fluid velocity) corresponding to the quiescent fluid of uniform 
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density po- 6 = is the definition of the normalization. We define the following hermitean 
scaled covariant derivatives Q: 

D k = Vd(-id yk + A k ) Vfc6{l,2} 

D = Vd(-idt + Ao) (2.2) 

Aq was introduced by Susskind Q to implement the constraint of the conservation of vorticity 
of the fluid into the action. By defining 

Di + iD 2 + Di - iD 2 , \ 

D= ^ , D> = * (2.3) 

and supposing that there is a constant magnetic field of strength B perpendicular to the fluid 
plane with vector potential 



A 



Bx 2 
2 

-Bx 1 
2 





(2.4) 



(one must not confuse A with (Ao, A\, A 2 ) as the first one describes the real magnetic field, 
B, and the second one the "position" of the fluid). The action (up to the second order in 
derivatives) for the fluid including the Maxwell term and the Chern-Simons term can be 
written as (2|: 

Smcs = J ^Tr{(-2[ J D , J D][ J Do^ t ] - [A [D, £>t]) + (2.5) 

+ 2A (~[D,Di] + l) D } (2.6) 

where 

■s 



is a rescaled time, g 2 = p0 . 



. _ eBO 1 ' 2 
A = 



m 

and e and m are respectively the charge and the mass of the electron. For simplicity we 
define: 

(2^V2) 
~ 2g 2 6 2 

The Chern-Simons term taken here corresponds to the one studied by Susskind M, and admits 
a smooth commutative limit. It is for this theory that in exact plane wave solutions were 
found and it was also found that the vacuum solution seems to be unstable. 

In order to reduce our system to a quantum Hall droplet we follow Polychronakos || 
and suppose x l ,x 2 and A^ \/fi £ {0, 1,2} to be N x N matrices where N is the number of 
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electrons. Like Polychronakos we subsequently introduce an iV dimensional vector \& which 
represents the boundary degrees of freedom. The boundary term introduced in the action by 
Polychronakos is — 2Y ^J^o^ which gives the total action: 



S = T J dt (Tr{(-2[D ,D}[D ,D^ - [D,D^][D,D^] 
+ 2A (-[£>,£>*] + l) A)} -2tftA)tf) 
Varying with respect to ^ we get the boundary constraint: 

D $> = 44> ii> = A fy 
Varying with respect to Dq we get the Gauss law: 

-[D, [D ,rt]] - [D\ [£>„,£>]] - K[D,rf) - 1) - tftft = o 
Finally, varying with respect to we obtain the Ampere law: 

[D Q , [D ,D]} + [D, [D, Z)t]] = x[D , D] 



Taking the trace of 2.10 gives: 



and therefore we can take \I/ = V NX I N — 1) where I N — 1) is a normalized vector. 



(2.7) 
(2.8) 

(2.9) 
(2.10) 
(2.11) 
(2.12) 



3. The Wigner crystal 



In this section we try to find a "vacuum" solution, that is a static solution (dt = 0) with 
An = (and therefore x l = y l ). The fact that dt is null imposes that | N — 1) is constant. 



Under those conditions the Gauss law ( 2.10 ) becomes: 

[£>,Dt] = t-N | JV- 1) {N- 1 | 
while the Ampere law ( [2.11| ) becomes: 

[D, | N - 1) (JV - 1 |] = D | JV - 1) (N - 1 | - | JV - 1) (JV - 1 | D = 



Replacing for D (with = 0) we see that equation 3J. is equivalent to: 



(3.1) 



(3.2) 



(3.3) 
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This implies in turn that projecting on the interior of the droplet (Vect{ | 0), | 1), | 2) 
N - 2) }) we have: 

[y\y 2 ] 



[x\x 2 ] 



(3.4) 



which endows the droplet with a noncommutative geometry (9|, [lOf| . 

A solution to the Gauss law (|3.1|) is the so-called Wigner crystal. That is: 

N-2 

D = a= Vn + 1 | n)(n + 1 | , (3.5) 



n=0 



where { | k) \ k G [| 0, N — 1 |] } is an orthonormal basis. The notation used here and 
throughout this paper, in general, u k G [| a, (3 |] " signifies that k is an integer that takes 
values between the integers a and (3 inclusively. The operators a and a) are very similar to the 
standard annihilation and creation operators with the notable exception that as \ N — 1) = 



(N — 1 | a = 0. Using 3.5 we have, projecting to the interior of the droplet 



\int. — ^ ~i~ 



„22 



\int. 



N~2 

(0/2)(Z>Z>t + rfD) \ int = (9/2) (2n + 1) | n)(n 

n=0 



(3.6) 



and so the electrons are spaced out in circles of radius oc \/n. This is the so-called Wigner 
crystal solution identified by Polychronakos M and it also corresponds to the "vacuum" 
solution of Susskind although he does not mention the connection explicitly. 

In our case of course it is not a solution as it is easy to check that it does not satisfy the 
Ampere law |3.2| . It was shown that this solution is unstable on the infinite plane Q but now 
we see that it is (classically) impossible on the droplet with our action. Actually, under our 
hypothesis (dt = and = ) there are no solutions, seeing as 3A and 3^ are incompatible. 



In fact, in the ordered base of 3.2 implies that 



D 



M 

if m 



(3.7) 



Where m G C, M is an arbitrary N — 1 x N — 1 matrix and is the null N — 1 vector. It can 



easily be seen that this form is incompatible with 3T . To find solutions we must therefore 
modify our hypotheses. 



4. Soliton solutions in the Hall droplet 



4.1 Arnaudon-Alexanian-Paranjape solution 



Following M we will look for static solutions dt = with the following ansatz: 
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N-l 



d = VeA = H n ) I n )( n 



n=0 

N-l 



D = ^ f( n ) \n)(n\a 

n =0 

$ = y/N\\N-l) (4.1) 



where a is the annihilation operator defined in 3.5 and 4>{n) and /(n) are to be determined 
by the equations of motion. The last term of the sum for D is superfluous as it vanishes. We 
note the following identities which hold for any function g(n), and we will use them in the 
analysis that follows: 

N-l N-l N-l 

M$3 9 ^> I n ^ n I)] = °( X^( n ) -g(n-l)] \ n)(n |) = (^2\g(n + l) - g(n)] \ n)(n \)a 

n=0 n=0 n=0 



N-l N-l N-l 

i ai dJ2 9( n ) I n )( n I)] = a H^2[9{n) -§{n + 1)] | n)(n |) = [g(n - 1) - g(n)] \ n)(n Da 1 " 

n=0 n=0 ri=0 

Notice that g{— 1) and ^(iV) are not defined, although they appear in the equations formally 
they are not in fact present because they are coefficients to states which are annihilated. 

Let us also calculate the "magnetic field", B, and "electric field", E, of the fluid: 

B = [D,D^ = [X)((n+l)|/(n) | 2 -n | f(n - 1) | 2 ) | n)(n \\ 

\n=0 J 

- (N - 1) | f(N-2) | 2 | N- 1)(N- 1 | (4.2) 
E = [D , D) = X /(n)Q(n) | n) (n | a (4.3) 

n=0 

where for n e [\ 0, TV— 2 |] Q(n) = 0(n) - 0(n + 1) and we define Q(N - 1) = <£(iV - 1). Let 
us also fix B{n) = (n \ B \ n). 



The equations of motion with the ansatz |4.1| gives for the boundary constraint [2.S| : 

= = 4 tf = v / iVA</>(iV - 1) => 0(iV - 1) = Q(iV - 1) = (4.4) 

This serves as the boundary condition for the first order difference equation relating 4>(n) 
to Q(n). The gauge field equations independently determine the values of the Q(nYs. The 



Ampere law 2.11 becomes 



N-l N-l 

^ f(n)(Q(n) 2 + B(n + 1) - B(n)) \ n)(n | a = A /(n)Q(n) | n)(n | o (4.5) 

n=0 n=0 
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while for the Gauss law 2.10 we have: 
N-l 



n=0 



2 £ ((n + 1) | /(«) | 2 Q(n) - n | /(n - 1) | 2 Q(n - 1)) | n)(n 



X[B - (1-N | N- 1) (N- 1 | ) 



(4.6) 



(4.7) 



Requiring that Vn G [|0, iV — 2|], /(n) / and replacing for B, 45 reduces to: 

Q(„)2 + ( n + 2) | f(n + 1) | 2 -2(n + 1) | /(n) | 2 +n | /(n - 1) | 2 = AQ(n) (4.8) 

Vn G [|0,JV - 3|] and 

Q{N - 2) 2 - 2(7V - 1) | f(N - 2) | 2 +(N - 2) | f(N - 3) | 2 = XQ(N - 2) (4.9) 

for n = N — 2. Then defining that /(iV — 1) = 0, since f(N — 1) does not actually appear in 
the definition of D we obtain one single defining equation: 



Q{nf + (n + 2) | f(n + 1) | 2 -2(n + 1) | /(n) | 2 +n | /(n - 1) | 2 = AQ(n) 
Vn G [|0, N — 2|]. We can rewrite this as: 

i(2Q(n)-A) 2 + V 2 ((n + l)|/(n)| 2 ) 



(4.10) 



A 2 
4 



(4.11) 



Where V 2 /i(n) = h(n + 1) — 2h(n) + h[n — 1) is the discrete one dimensional Laplacian. 



Furthermore one easily shows H by induction that 4.7 reduces to: 



(2Q(n) - A) | f{n) | 2 = -A 



(4.12) 



Vn G [|0, N - 2\]. Then replacing |4.12| in [4.11| we get: 

2 



1 



A 



+ V 2 ((n + l)|/(n)| 2 ) 



A~ 
4 

.2 _ A 2 



(4.13) 



Vn G [|0, N — 2|]. Finally, defining u n = (n + 1) | f(n) \ and g = ^- we get the following 
recursion relation: 



u n+ i - 2u n + n„_i + g 



2 (n + l 



9 



(4.14) 



With the boundary conditions that n_i = (-1+1) | /(-l) | 2 = Oandujv-i = N \ f{N-l) | 2 = 
0. Which is exactly verifies the suggestion in M. 
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4.2 New soliton solutions 



Our choice of ^ oc| N — 1) is quite arbitrary, and one might wonder what would have been 
different had we chosen differently. This is what we will now investigate and we will see 
that it leads to new soliton solutions. Take 



N-l 



n) 



n=0 



We know from |2.12 that 



N-l 



El A, 



(4.15) 



n=0 



Keeping the same ansatz (|4.1|), all our equations remain unchanged until 4/7 which becomes: 
N-l 

2 ^((n + l)| f(n) | 2 Q(n) - n | /(n - 1) | 2 Q(n - 1)) | n)(n | 

n=0 

JV-1 JV-1 

= XB - XI + ^ AmA ™ I m )( n I 



(4.16) 



m=0 n=0 



Since both the left-hand side of this equation and B are diagonal we must have m ^ n 
X m X* n = 0, which in turn implies that * oc| M) with M £ [\0,N - 1|]. Thus we take 
\P = \J NX | M) for some fixed but arbitrary value of M. Then the Gauss law ( f4.16| ) becomes: 



N-l 

2 ^2 (in + 1) | f(n) | 2 Q(n) -n \ f(n-l) | 2 Q(n - 1)) | n)(n 

n=0 

= XB - Al + iVA | M)(M | 



(4.17) 



By ^1] we now have <j){M) = and </>(iV — 1) is no longer zero. This is now the modified 
boundary condition that is used in the difference equation relating the ^(ra)'s to the Q(n) 7 s. 
One easily shows by induction that equation 4.17 reduces to : 

(2Q(n) - A) | f(n) | 2 = 



-A 

NX 
n+1 



VnG [|0,M-1| 



A VnG M,JV-1 



(4.18) 



however consistently still imposing f(N — 1) = 0. Equation 4.11 is the Ampere law and 
remains unchanged. Inserting 4.18 in [4.11 yields: 



4 



A 



4 VLfWF 

1 / (N- 
4 Un+l)|/(n)| 



+ V 2 (Jn+l) | f(n) | 
^^-^Vv 2 ((n+l)|/(n 



Vne [|0,M-1|] 
VnG [|M,iV-2| 



(4.19) 
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And so we get the modified recursion relations: 

( n _|_ i\ 2 

u n+ i - 2u n + u n _i + g 2 = g 2 Vn E [|0, M — 1|], 

\ Un J 

(N — 1 — n\ 2 
= fl 2 VnG [\M, N - 2|], (4.20) 

with the same boundary conditions as before, U-i = ujy-i = 0. Thus we have in all N 



solutions for the ansatz 4.1. We can prove that \/g 6 M and VM £ [|0, iV — 1|], there exists 
a unique no > such that U-i = = itjv-i and u n > OVn 6 [|1, iV — 2|]. This is intuitively 
obvious and we relegate a rigorous proof to the appendix. Furthermore, due to the obvious 



symmetry of 4.20) and the symmetric condition u_i = mat_i = we have that: 

UM=m,n=k = u M=(N-l-m),n=(N-2-k) (4-21) 



One might be lead to think (because of 4.21) that the solutions M = m and M = N — 1 — m 



are the same solution with only a permutation of the basis vectors. This is in fact not the 



case because this permutation (or gauge transformation) does not leave a (as defined in 3.5 ) 
invariant. The energies of these solutions, however, turns out to be degenerate. Nevertheless, 
we have N distinct solutions. The recurrence relations 4.20] cannot be solved analytically, in 



figures 1-4 we show some numerical aspects of the solutions. 
Let remark also that by fixing 

\& = NXe^ | M) (4.22) 



instead, we can shift the <^>(n)'s by an arbitrary constant a. The boundary condition for 
the </>(n)'s changes correspondingly. We can also translate our solutions to the gauge Aq = 
albeit time dependent solutions. Our static solutions relied only on the following commutation 
relations: 



N-l 



[D ,D] = d{n)Q{n) \ n)(n | a 



n=0 
N-l 

[D ,D^] = ^ -d(n)*Q(n)a j \n)(n\ (4.23) 

n=0 



where 



N-l 

D = \ d{n) | n)(n \ a. 



n=0 



The second commutation relation is a consequence of the first one if the Q(n)'s are real (which 
they are, according to 4.18[ ). But these commutation relations can be obtained with a time 



dependent D with Aq = (and thus Dq = —iyOdt) by taking d(n) = f(n)e ^ . The other 
equations are obviously fulfilled since they depend only on the modulus of d(n) (or f(n)). 
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Figure 1: Behavior of the solutions of the re- 
cursion equations 4.20 for the different values 
of M. The plot shows u n for N = 101. 
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Figure 2: Behavior of the solutions of the re- 
cursion equations 4.20 for the different values 
of g. The plot shows u n for TV = 101. 




Figure 3: Behavior of uo as a function of M. 
Numerical solutions for g = 1 and N = 101. 
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Figure 4: Zoom of Figure 3 around the value 
of M = 50. 



5. Properties of the solitons 



5.1 Energy 



The Hamiltonian corresponding to action 2.8 is: 

H = ~ffr[ - 2[D ,rt}[D ,D] + [D, D*][D, rf] + D } (5.1) 



The term TtDq = ^2m=o ( t ) { n ) vanishes in the time dependent case and can be put to zero 



in the static case by adding a phase to according to |4.22 , So, we will henceforth suppose 
TtDq = 0. The kinetic energy is. 



T = -^=Tr{ -2[D , rf][Do,D]} = l-^Tr&E 



(5.2) 
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The potential energy is 

nun ^-y-" 

V = -j=Tr{[D,rt][D,Di]} = -jfr& 
For simplicity define H = = then replacing for D and in 5.2 we get 

JV-2 N-2 

T = 2~ £ (n + 1) | /(n) | 2 Q(n) 2 = 2~ £ Q(» 



(5.3) 



[n) 2 u n . 



(5.4) 



n=0 

In terms of the u n 's this means: 

M-l 



ra=0 



2 2 2 



(5.5) 



n=0 " n=M 

Using [O^ we can rewrite this last equation in a form more convenient for numerical analysis: 

N-2 



T = 2H { 9 2 u n ~ u n V 2 u n ] +g 2 (M 2 + (N-M- l) 2 ) 



(5.6) 



n=0 



Equation 4.21 tells us that the sum in 5.6 is the same for M = m and M = N — 1 — m and 
it is easy to verify that the additional constant also has this symmetry. Therefore 



T(M = m) = T(M = N-l-m). 



(5.7) 



Nevertheless, the solutions for m = M and for M = N — 1 — m are not gauge transforms of 
one another. As is clear, the permutation symmetry does not commute with the operator D, 
and the distribution of the "magnetic" field is different. 



Likewise, by replacing for B in |5.3| , we get the following expression for the potential 
energy: 

N-l 

\2 

n=0 



V = E ^2(u n - u n -iY 



(5.8) 



Due to 4.21 we have that: 



V(M = m) = V(M = N-l-m) 



(5.9) 



And so, by adding |5.6| and |5.8| and taking into consideration that u-± = itjv-i = we get the 
following expression for the total energy: 



N-2 



N-2 



E = E \j2 u n^ 2 u n + g 2 {4 un) +n({N-1)-M 



(5.10) 



n=0 \ n=0 

In figures 5, 6, 7 we show some numerical aspects of the energy for various values of M. 
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"Kinetic Energy" 




"Potential Energy" 



Figure 5: Behavior of T as a function of M 
in units of S. Numerical solution for g = 1 
and N = 101. 
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Figure 6: Behavior of V as a function of M 
in units of S. Numerical solution for g = 1 
and TV = 101. 
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Figure 7: The energy of the solutions for different values of M. Numerical values for g = 1 and 
N = 101. 

5.2 "Magentic" field and flux 

The total "magnetic" flux for any of the solitons is strictly zero. This is a simple consequence 
of the fact that B = [D,D^\ hence the total flux 



total 



TrB = Tr[D,D^} = 0. 



(5.11) 
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However, if we compute the flux up to a state | R — 1) which corresponds to the flux localized 
in a radius ~ \/R — 1, we find 



R-l R-l 



&R-1 



(m I i? I m) = (m | [D, D^] \ m) 



(5.12) 



m=0 m=0 



(u m - U m -l) = u(R - 1) = Rf(R - l) 2 



(5.13) 



m=0 



Comparing with Figure 1, we see that for different values of M, the flux is concentrated over 
different regions. For small M, there is a positive flux tube at the origin, surrounded by 
a wide negative flux annular region, while for large M there is a wide cylindrical region of 
positive flux surrounded by a localized annular region of negative flux near the boundary. 
The total flux always vanishes since u(N — 1) = 0. 

6. Conclusions 

In this paper we have found soliton solutions of the noncommutative Maxwell-Chern-Simons 
theory, a modified version of the pure noncommutative Chern-Simons theory studied by 
Susskind Q and which he proposed as the theory of the quantum Hall effect. In a previous 
work ^ we had studied the modified theory on the infinite noncommutative plane, and we 
had found plane wave solutions. Vortex-like soliton solutions were however more elusive and 
it seems that, in fact, the Susskind vacuum is unstable to the perturbation by the Maxwell 
term. In this paper we restrict the Maxwell-Chern-Simons theory to a finite droplet adding 
boundary degrees of freedom, as proposed by Polychronakos HJ. We find non-trivial soliton 
like solutions to the equations of motion. It appears that the energy is minimized for the most 
symmetric solution which corresponds to a positive "magnetic" flux cylindrical core over the 
center half of the droplet surrounded by a negative flux region over the outside half. It should 
be noted that this "magentic" field is an analog of the usual Maxwellian magnetic field and 
in fact here it actually corresponds to vorticity in the fluid velocity. The total flux is always 
zero. 

The model considered is directly related to the model studied by Polychronakos || and by 
Hellerman and Van RaamsdonkQ, which is shown to be equivalent to the Laughlin theory [j?]] 
for the quantum Hall effect. Their model however, contained a harmonic oscillator potential 
the sole aim of which was to break the translational invariance of the noncommutative plane. 
We do not consider this harmonic potential. Instead we add the next order corrections, after 
the Chern-Simons term, in the energy of the velocity of the quantum fluid, the Maxwell 
term. We find that adding the Maxwell term a more satisfactory modification of the pure 
Chern-Simons theory, at least the Maxwell term does measure the actual energy of the flow in 
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the fluid. With the Maxwell term added we find that there exist multiple classical solutions 
which are essentially vortices located at the origin superposed with a background of opposite 
vorticity. Every physical sample of quantum Hall material exhibits a transition to the Hall 
insulator at high enough magnetic field ||. This transition exhibits a perfect duality in the 
current to voltage curves just above and just below the transition. The current flow just above 
the transition is attributed to charge transport by vortices while just below the transition it 
is attributed to particles. The vortices that we have found should model these very physical 
vortices of the quantum Hall system. 

There are several avenues for future study. First one should find the analog of plane wave 
excitations in the model. Then as a function of the parameters of the model (the external 
Maxwellian magnetic field, the density (related to the noncommutativity parameter 9 and the 
coupling constant g (related to the relative strength of the Maxwell term to the Chern-Simons 
term) it would be very interestiing to find a critical theory where the plane waves and the 
vortex-like solitons become degenerate. At this point there should be a phase transition in 
the behavior of the theory. This could afford a description of transitions between plateaux 
or the transition to the Hall insulator that are observed in any experimental Hall system. 
Finally, from a more mathematical point of view, it would be interesting to know the full 
modulli space of solutions to the equations of motion that we have studied. The stability and 
other properties of the solutions should be classified. 
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A. Proof of the existence and the uniqueness of uq for each M 



In the following analysis we will assume that u_i = and that the u n 's simply satisfy the 
recurrence relations 4.20) . Then we will prove that there exists a unique value of uq for which 
the recurrence relations 4.20 and the boundary condition itjv-i = are satisfied. 



Definition Let U = { x G R + | u = x u n > 0, Vn G [|0, N - 1|]} 
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Definition Let = { x G R+| n = x => u n > Vn G [|0, TV - 2|], and n^v-i > 0} 



Lemma A.l 0^0, [iV, oo) C O. 



Proof Take x > N. We will show that x G fi. We use induction to show that no 



u,. 



> N, Vm G [|0, AT" — 1|], and u m > u m -\. First, if uq = x, then evidently uq > N and 
u o > = u_i- Thus the induction hypothesis is valid for n_i and no- Next we assume that 
the induction hypothesis is valid for each integer k less than or equal to a fixed n G [|0, N — 2|]. 
That is, VA; G [|0, n|] we assume Uk > nfc_i > A 7 '. With this assumption we will prove that 
Un+i > u n > N. By 4.20 we have 



U n +1 



-1+^1 



X(M, n) 



(A.l) 



J(M, n) 



where 

(n+1) 2 for n G [|0, M — 1|] 

(N-(n + l)) 2 for n G [\M, N — 2\] 
Evidently, < X(M, n) < N. By the induction hypothesis we have u n > n ra _i > N. Hence 
H^Bi < i. Thus, using the recurrence relation A.l , we obtain n n +i > u n > N > 0. 
Therefore, by the principle of mathematical induction, u m > N > 0, Vm G [jO, JV — 1|], if 
u G [AT, oo). Therefore [N, oo) C fi ^ 0. □ 



The n n 's are smooth functions of uq on f2, Vn G [|0, AT — 1|], since they never vanish on 
ft. They are in fact simple rational polynomial functions of Uq. Evidently, the only way to 
introduce a singularity is through the recurrence relations 4.20[ , which become singular only 



when any of the n n 's vanish. This is of course also true on f2, since only njv-i may vanish on 

d U n 
d Uq 



f2, which only introduces a singularity in un, which is irrelevant. Thus 4-^ is well defined 



onfiVnG [|0,JV-1|]. 

Lemma A.2 Vn G [10, N — 111, 4-^ > and > d -^^ on fi. 

u > u' a uq a uq a uq 

Proof We will again use proof by induction. First for n_i and no, evidently = 1 and 
d , u ~ 1 = since n_i =0. Hence 4-^ > d , u ~ 1 ■ Next we assume that the lemma is true 

a uq 1 d uq d uq 

Vfe G [|0, n — 1|] for some fixed n G [|0, A" — 1|] and then prove we that it is true for n. This is 
simply obtained from the recurrence relation 

n— 1 



d u n d n n _i d« n _2 2 I(M,Ti-l)rf ?t. 



_2— —--r — + 2g , 

a no a no a no n n-i " u o 

> , + , r^) > , > (A.2) 

a no a no a no a no 
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using the induction hypothesis and that X(M, n — 1) > 0. Thus 



dUn > dUn-l >Qyne [|Q) ^ _ n (A3) 

a uq a uq 



Definition Let /i = inf(f2). 



is a strictly positive real number since for example, f2 C (<5, oo) where (5 is defined by the 
value of uq = 5 > which renders u± = 0, that is 

„, = „ = « + ," (l-3£°>). (A.4) 

Using our two lemmas, if we start with uq > N we know we are in fi, and if we now reduce 
the value of uq to 5, we know that the value of u\ will decrease monotonically to zero (we will 
in fact encounter singularities in the other n n 's already, but for the present purposes these do 
not matter), hence we are no longer in 0. Thus 0, C (5, oo) and [i > b > 0. 
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Proposition A. 3 u is the unique element o/K such that uq = \i u n > 0,Vn G [|l,iV — 2|], 
and un-i = 0. 



Proof First of all we see that n = inf(O) ^ Vt. This is because all of the u n 's are smooth, 
continuous functions of uq on f2. Thus if n were in Q, then for tig = n, 3 e > 3 u n > 
e Vn G [|0, iV — 1|]. Hence by continuity of the functions u n there exists a neighbourhood of 
[i for which u n > which contradicts the hypothesis that n = inf(fi). This in turn tells us 
that for uq = n we have one of the three following possible cases: 

1. 3n G [|0, N - 1|] such that u n = 

2. 3n G [|0, N - 1|] such that n n < 

3. 3n G [|0, iV — 1|] such that n n is ill defined. 

The third case is clearly a subcase of the first, for if k is the smallest number such that 
is ill defined, k > 2 since uq = \i ^ and hence u\ is not ill defined, so the first possible ill 
defined u^ can be U2- However, by assumption u^-i and u^2 are well defined and then by 



the recurrence relation 4.2C , 



Uk = 2u fc _i - u k _ 2 + g i 



J(M, fc - 1) 
u l-i 



But this is ill defined only if n&_i = 0, thus we are necessarily in the first case. 

The second case can also be easily ruled out. By continuity of the u n 's, if u n (fj,) < for 
a sufficiently small e > we still have u n (/j, + e) < but fi + e£Qas/j> = inf(O). This is in 
contradiction with the definition of O. 

Therefore only the first case is possible. Suppose now that for some m G [|0, N — 2|] and 



u m (uo = fi) = 0. Then by the recurrence relations [4.2C 



lim n m+ i(n ) = -oo. 

Then for a finite e > 0, we have fj, + e G O,, but u m +i(n + e) < which is a contradiction. 
Therefore n m (n = |i)/0VmG [|0, N — 2|]. 

We are therefore forced to conclude that when no = /i =^ itjv-i = and n n > 0, Vn G 
Nl,iV - 211. 



Moreover lemma |A.2j implies that njv-i is an injective function of no on Q,, proving the 
uniqueness. Therefore \i is the only solution. □ 
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